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1 Introduction 



Any physical system can be studied by using different sets of variables which 
yield, of course, mathematically different descriptions. Different choices of 
variables must however be physically equivalent, in the sense that they should 
describe one and the same system. An extreme manifestation of this property 
appears in models which can be defined in terms of either fermionic or bosonic 
variables. The equivalence between these two formulations is made explicit by 
the so called 'bosonization rules', that map fermionic into bosonic variables. 

Bosonization turns out to be a very useful tool indeed in order to under- 
stand and in some cases even to solve some non-trivial interacting Quantum 
Field Theory models in two spacetime dimensions. It is an interesting fact 
that there is no theoretical obstacle to the extension of this procedure to 
higher dimensions. Indeed, there has recently been some progress in the 
application, although in an approximated form, of a path integral bosoniza- 
tion procedure to theories in more than two spacetime dimensions 
dealing with both the Abelian and the non-Abelian cases. An important 
difference between these works and two dimensional bosonization |[T2|-|T^ is 
that, in the latter case, the fermionic determinant can be evaluated by per- 
forming a decoupling transformation of the fermions |T^]-[^. The redefined 
fermionic fields (in the massless case) are free, and the effects of the inter- 
action manifest themselves through the existence of an anomalous Jacobian 
which depends on the gauge field. We remark that no such a thing had yet 
been suggested for the higher dimensional case; the present work deals with 
such an extension. In this paper we shall complement the usual path integral 
bosonization technique by showing how to redefine the fermionic fields in the 
path integral, in order to decouple them from the gauge field. This step, 
we believe, fills a gap in this approach to bosonization in higher dimensions, 
rendering the whole process entirely analogous to its well known two dimen- 
sional counterpart. The main point of our work is that it gives, contrarily 
to previous bosonizations in 3 dimensions, the bosonization with an explicit 
expression for the fermionic fields in terms of bosonic ones and free fermionic 
fields. This procedure is moreover generalizable to higher dimensions. 

Let us briefly review, for the Abelian case in 3 dimensions, the usual 
approach to bosonization, which starts from the definition of the generating 
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functional 2 {J), 

Z{J) = j V^V^ exp - J d^x^{^+ / + M)* 

= det [^+ / + M] = exp [-W{J)] (1) 

where W{J) is the generating functional of Euclidean connected correlation 
functions of fermionic currents. The bosonized action 5'(,os(^) is then given 
by a generalised functional Fourier transform of Z{b) = exp[—W{b)]: 

eM-Sbos{A)] ^ j Vb^e-x^ -W{h)+i j d'xe^.xbi.d.Ax , (2) 

and the bosonization rule for the fermionic bilinear ^7^^^ is 

f'(x)7^^(.T) ^ iei,,xd.Ax{x) . (3) 

The fact that the functional W{J) proceeds from fermionic matter field 
seems to be here a matter of no relevance, the crucial property invoked being 
gauge invariance of the fermionic determinant instead. This is quite different 
to what happens in two dimensional path integral bosonization, where the 
fermionic fields are redefined in terms of new, decoupled fermions, through 
transformation functions that depend on the bosonic gauge fields. In this 
paper, we give explicit transformation functions that achieve the same goal 
as in the two dimensional case. 

A difficulty renders the bosonization procedure in more than two space- 
time dimensions non-exact. It is our inability to compute exactly a fermionic 
determinant in those cases. This allows us to determine only approximated 
bosonized actions, which usually result from some perturbative or low-momen- 
tum expansion. We shall see that the same can be said about the decoupling 
transformations, one can consistently decouple the fermions up to some order 
in the relevant coupling constant of the model being considered. 

The structure of this paper is as follows: In section 2 we explain the 
mechanism of decoupling on a general footing, discussing both Abclian and 
non-Abelian cases, in the massless and massive situations. The application 
of the general results to the path-integral bosonization of the three dimen- 
sional Thirring model is studied in section 3, and in section 4 we explain the 
differences and similarities between our approach, when restricted to the two 
dimensional case, and the one which makes use of the anomalous (Fujikawa) 
Jacobian method. 
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2 Decoupling the fermions by a change of 
variables 



2.1 Massless case 

The generating functional Z{A) for massless fermionic fields in the presence 
of an external gauge field A in a D-dimensional Minkowski spacetime is 
defined by 

Z{A) = yp^rD^e'^^^*'*'^) (4) 

where 

JP = -i^D^ D^ = + eA^ (5) 

and 

_ _ f ^-^M in the Abelian case , ^ 

~ 1 A'^^Ta in the non Abelian case ' 

where Tq are (anti-hermitian) generators of the Lie algebra of the non- Abelian 
gauge group, and both A^ and are real. Dirac's 7-matrices are assumed 
to satisfy 

{7m > 7 J = 2c/^^ 7o = 7o l] = -Ij ■ (7) 

We shall show that, by a change of variables, one can decouple the 
fermions from the external field A. Then we redefine the fermionic fields 
^{x) in terms of new ones x{^)) follows 

^(x) = ^(er'4)x(a;) 

^>{x) = x{^)T{er'f) (8) 

where 

= 7o[^(er'4)]So (9) 

and JF is a function which, we shall assume, is defined by a power series 
expansion in e 

00 

= 1 + E ""e"(r'4r, (lo) 

n=l 



^We shall later on explain the analogous procedure for the Euclidean spacetime case. 
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with real coefficients The fact that we shall understand the transfor- 
mation (P) as given by the series ([T0|) is the reason why we call the whole 
technique 'perturbative'. 

It is important to note that the operator introduced in (P) must be 
understood as acting on both ^ and the fermionic field Xi iiot on / alone. 
Namely, 

^(er'4)7^^(e(r'4)), (11) 

where the outer parenthesis denotes functional dependence, while the inner 
one delimites the action of the operator . This fact and the reality of the 
a„'s implies that the transformation formula for ^>{x) may also be written 
more explicitly as 

^{x) = x(x)^(e4r') (12) 

and this property turns out to be crucial in the decoupling. We then see that 
the new fields will have a free action as long as the equation 

j d''xxHe4r')^V>Her'^)x{x) = jd^'xxi^xix) (13) 

may be satisfied. This implies 

^(e4^-^) =^ (14) 

which, by using the power series expansion (p!0D , yields a set of relations for 
the unknown coefficients a„. They may be written as 

oo 

= e(l + 2ai)/+5]2K + a„+i)e"+ir(") + 

n=l 

oo n oo n—1 

Y^iY. «n+l-m«m)e"+^r(") + ^(E «n-„^«„^) e^+^T^") (15) 
n=l m=l n=2 m=l 

where 

t(") =4{^-^^y = {j^^-^Y =4^-^ ^-^ ^ (16) 

namely, the operator T'-"-' consists of a product of alternating factors and 
^, starting and ending with an ^, and containing n factors altogether. 
The crucial property that makes the decoupling possible is that any redefi- 
nition like (H), (even with arbitrary coefficients a„), modifies the action by 
adding a series of terms, which are always proportional to one of the T'^'^^'s. 
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Equation (|T5|) implies the recurrence relations 



3 1 """^ 

ftn+i = -T"n - -J^i'^k+i + ak)an-k Vn > 1 (17) 

k=l 



plus the initial conditions 
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«i = -- «2 = - . (18) 



The solution to the recurrence relations (|T7D, with conditions ([T8|), is 



Or 



_ 1 

2 
T2, 



where 



H)H-i)---H-^ 



(19) 



(20) 



Thus we conclude that J-'{x) can be thought of as the power series expansion 
of the function (1 + z)^2 around x = 0, with a radius of convergence equal 
to 1, 

H^) = E f ~M = (1 +^)"^ ' 1^1 < 1 ■ (21) 

n=0 V 'W 

Condition |x| < 1 is equivalent to saying that we are in the perturbative 
regime. Whenever we write, henceforward, jF(x) = (l + a;)~2, we shall have 
in mind its power series expansion as given by (^TJ). 
When performing the transformations 

^{x) = {l + e^-^ 4)-^ x{x) 

^{x) = x{x){l + e ^r'rK (22) 

we can then affirm that the new fermionic action shall be free, but, on the 
other hand, there will appear a Jacobian J{A), due to the change of (Grass- 
mann) variables, 

Z{A) = yp^D^ e^^^(*'*'^) 

VxVxJiA)e'^^^^''''^'' (23) 
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where, 

Sf{x,X;0) = j d''xx{x)i^x{x) (24) 

and 

J{A) = JJA)J^(A), 

UA) = {det[(l + er'4)"^]}"' 

J,(A) = {det[(l + e4r')-^]}"' • (25) 
We remark that these Jacobians are different from one, but nevertheless not 
anomalous, in the sense that, as the field transformations are not unitary, 
they shouldn't be equal to one. 

From the two self-evident relations 

det^ = det ^det(l + e^~^4) 

det ^ = det(l + e 4^"^) det ^ (26) 

we deduce 



Recalling (0), we see that, as it should be, multiplying the Jacobian by the 
free determinant we reconstruct the determinant of p. 

The use of transformations (^) would, in principle, allow us to decouple 
the fermions completely from the external field A. This, however, is not 
possible to use in the apphcations, since it would require the evaluation of the 
Jacobian (|27| ) exactly, something not possible in more than 1 + 1 dimensions. 
We may, however, still use a finite order approximation to (p2D, in order to 
achieve partial decoupling. For example, we can perform a transformation 
that corresponds to the truncation to second order in e of (|2^). This will 



lead to an action for the new fermions containing interaction terms of order 
three or higher; namely, 

^{x) = (l - ^er' 4 + le'r' 4) Xix) 

^{x) = x{x) (1 - le + 4r') (28) 



and 



j d''xx{i^+0{e'))x. (29) 
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The use of the decouphng transformations in the Euclidean spacetime 
case could of course be implemented by first rotating to Minkowski space- 
time, then performing the transformations and at the end coming back 
to Euclidean. The new fermions will of course remain decoupled after the 
rotation. If one wanted to implement the transformations directly in Eu- 
clidean spacetime, one should face the problem that the Dirac matrices have 
different hermiticity properties than in Minkowski spacetime, and that spoils 
the transformation rule for the adjoint spinor, since under the transformation 

^(x) = T{er^ ^)x[x) (30) 

it is still true that 

^{x) = x{x)T{er^f) (31) 

but 

7o 



^(er'4)lSo = T{el~f^) ^ HeW') (32) 



where, for any Euclidean vector t>^, we denote = (— Vi, — f2, fs). 

Although not necessary in principle, one can, however, still give an explicit 
decoupling transformation in the Euclidean case. To that end, we employ 
the trick used in where it was introduced in a slightly different context. 
That trick amounts to make in the Euclidean path integral representing the 
fermionic determinant 

Z{A) = yp^ri)^ e-/'^°"*(")^*(") = det^, (33) 

the change ^{x) — ^ 'I'(a;)73, which of course does not change the value of the 
determinant 0. Then one see that the fermionic action is 

Sf = J d^x^{x)^3JZ>'^{x) (34) 

and, under the decoupling transformations 

^(x) = {l + e^'^ 4)-'^ xix) 

^{x) = x{x)-/s{l + e ^r'rhs , (35) 



^One could of course have just said that one wanted to evaluate the determinant of 

73^. 
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the new action is 

Sf ^ J d^xxix)-f3 ^x(x) (36) 

as can be easily checked. Note that the presence of 73 in the new action is 
as irrelevant as in the old one. To avoid writing the 73, when we use the 
Euclidean formalism in section 3, we work with "^"^ instead of ^, 

Sf = j (Px^\x) ]p^{x) (37) 

and give the transformations in terms of those fields. The procedure is, of 
couse, equivalent. 



2.2 Massive case 

The massive case admits, within our decoupling approach, two different treat- 
ments. Beginning from the massive fermionic action 

Sf{^,-^;A) = J d^x^{x){ip - M)-^{x) , (38) 



one possibility is to perform the analogous of the redefinitions (22) but in 



eluding the mass term into the inverse of the free Dirac operator, i.e., 

^(x) = [1 + + iM)'^ xix) 

^(x) = x{x) [I + e + iM)-^]-^ . (39) 

This way of transforming the fermionic fields leads, after some algebra that 
is essentially equal to the one of the massless case, to the decoupled fermionic 
fields X cLnd x, with a massive free action 

Sf{x, x) = I d^'xxix) {i^- M) x{x) , (40) 

and to the Jacobian 

_ det(^ + ^M) _ 
-^^^^ - det(^ + zM) ' 

this way of decoupling, applied in its second order approximated version, 
is the one we use in the next section to study the Thirring model in 2 + 1 
dimensions. 
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There is, however, another way of transforming the fermions, which is 
in fact more similar to the one used in two dimensions. It consists in using 
the transformations ([2^), corresponding to the massless case, for the massive 
action. The new fermions are no longer decoupled, but all the interactions 
are proportional to the mass, and the new action is a suitable starting point 
for performing an expansion in powers of the mass, as it happens to be the 
case in two dimensions 

-M j d'^xx{x)[[l + e{^ + iMY^ 4][l + e4{^+iMY^]Y'' X{x) ■ (42) 

In Euclidean spacetime, as for the massless case, the continuation to 
Minkowski prior to decoupling is a possible path. However, one can also 
work directly in Euclidean spacetime as explained for the massless case. One 
then uses the 73 trick. If the first version of decoupling (^) is used, due to 
the presence of the mass term, we must also make a "smooth continuation" 
(in the sense of [^) to imaginary masses before decoupling, and then rotate 
back to real values. In perturbation theory, there is no obstruction to that. 
Then the transformations in terms of \E' and \E'^ read 

^{x) = [1 + + M)-^ 4]-^ x{x) 
^t(2;) = xHx) [I + e + M)-^]-^ . (43) 

where M is regarded as imaginary. 

We end up this section by making a consistency check, valid for both the 
massive and the massless cases. Namely, that the expectation value of the 
fermionic current ^'7^\E', when the fermions are written in terms of the x 
and X) yields the proper result. For the sake of concreteness, we deal the 
case of the massive action (^). The expectation value of that current, when 
evaluated in terms of the decoupled fermions is 

{^{x)^^^{x)) = 

{X{x)[l + e4{^ + 2M)-i]-57^[l + e(^ + tM)-' 4]-^ ^(x)) , (44) 

or 

{^ix)-f^^ix)) = 
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-Tr j d^yd^z{{y\[l + e4{^ + iM)-^]-'2\x)j^ 



(45) 



where the usual notation (x| ■ ■ ■ |?/) denotes kernels of functional operators in 
coordinate space. 

Since the new fermions are free, 



ixizMy)) = {z\{^ + tM)-'\y) . 
Using this, and after some straightforward algebra, we obtain 



(^'(x)7^^(x)) 



-Tr 



'^,,{x\i^+ie4 + iM)~^\x) 



(46) 



(47) 



which is the proper expression for the expectation value of the current in an 
external field. Also, (^) can of course be written in terms of the fermionic 
determinant: 



(^(a;)7^^(x)) 



eSA^{x) 



lndet(z p - M) . 



(48) 
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3 The Thirring model in 2 + 1 dimensions 



We shall here follow the same strategy as in for the 1 + 1 dimensional 
Thirring model, but with the necessary changes due to the fact that we are 
now in 2 + 1 dimensions, and the decoupling is, in consequence, non-exact. 

We define the generating functional of current correlation functions for 
the Thirring model in 2 + 1 Euclidean dimensions by 

Z{J) = Jv^^V^ expS^- J d^x[^\^+ ^+ M)^ + ^g{^^-f^^f]^ (49) 

where g is the coupling constant and is a source introduced in order 
to generate fermionic current correlation functions. In order to render the 
fermionic action quadratic, we introduce an auxiliary vector field A^, such 
that the generating functional now reads 

Z{J) = J V^^V^ VAf, 

xexpi^- J d^x[^\^ + t^4+^+M)^ + ^{A^y]^ . (50) 
Making now a shift in the vector field A^ ^ A'^ = i^A^^ + J^, 

Z(J) = y"r'^tx?^r'A;exp|-y"d=^x[^^(^+/ + M)^ - ^(A'- J)2]| . 

(51) 

Now we proceed to make a decomposition of the field A into the curl of a 
vector field plus the gradient of a scalar field c/?, 

= i{e^uxdu^x + d^'^) . (52) 

We note that we can shift the vector field $^ by a gradient without affecting 
the configuration of A!^. This freedom is due to the fact that A!^ has three 
components, while on the right hand side of (0) there appear to be four 
(three from <l>^ plus one from (/?). This apparent contradiction is solved by 
the above mentioned 'gauge invariance' under transformations of which 
allows we to impose a gauge fixing condition on in order to leave only 
two components. We choose the Lorentz condition 

d-^ = 0, (53) 
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which shall, of course, be included in the generating functional. We then get 
an expression for Z{J) as an integral over the fermions and the bosonic 

fields (/9 

Z{J) = Jv¥V^V^i,Vip6{d-^) 
exp |- J d^x[^^^ + h^^e^^xdu^x + i + M)^ 

-^{ief,uxdu^x + idf,ip - J^f]^ . (54) 
By a redefinition of the fermions, 

qj\x) = ^'\x) e'^(") (55) 
we entirely decouple them from the scalar field Lp: 

Z{J) = J P^'^ V^' V<^^ Vip d{d ■ $) 
exp |- 1 d'x[^'\^ + i^^e^uxd.^x + M)^' 

(56) 

At this point we make use of the results of the previous section, to per- 
form a new redefinition of the fermions, corresponding to (0), this time to 
decouple them from $^ 

m'{x) = [l + {^ + M)-\^^)]-'2x{x) 
^'\x) = x\x)[l + {^<p){^ + My']--2 , (57) 

where we have used 

tlf^e^uxdu^x = (^f) (58) 

which holds as a consequence of the Dirac algebra in 3 dimensions when 
the Lorentz condition for $ is used. Transformation (0) introduces a $ 
dependent Jacobian, 

Z{J) = J Vx^VxV^f,VipJ{^)5{d-^) 
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expj- J d^x[xH^+M)x + ^^^(-a^V + ^M-^.)*. 

+ - Y/' + -gM9,^ + eM.Aa.$A)]| (59) 

where 

J($) = det[l + + M)-\^ ^)] . (60) 

As already explained, it is not possible to evaluate the Jacobian J($) exactly. 
We shall instead use the quadratic approximation of 0, which yields for this 
object 

J($) = exp[-Ty($)] (61) 

with 

W^($) = I d^x^^[C^V+ + C-V-]^,^, (62) 

where the C± are scalar operator functions and the V± are two of the 
three orthogonal projectors used in the calculation: 

\ d d i d d 

[P±]^u = -{S^u - ± -^==e^xudx) Qt^v = • (63) 



We shall also need the transverse projector 

--Q2 ■ 



r = V+ + V. = S,.-^ . (64) 



Integrating out the field tp, 

Z{J) = J Vx'VxV<^^6{d-<!>)expS^- J £x[x\^ + M)x 

+ W{^) - - ^J^.V^.Ju + ^Jf^e^uxdu^x]^ . (65) 

The term quadratic in the source can be, as in 1 + 1 dimensions, linearized 
by the introduction of an auxiliary field. In 2 + 1 dimensions it has to be a 
vector field and the new expression for Z becomes 



Z(J) = J Vx^VxV^^Vi^5{d-^) 
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exp |- 1 d^x[xK^ + M)x + W{^) - + ^{d ■ 0' 

+ ^^,d'V,.C, + tJ,e,,xd,{h^ - PCx)]^ (66) 

where the constants a and (3 satisfy the relation: 

/?2 1 

- = - • (67) 

a g 

and A 7^ is arbitrary. 

We then introduce two vector fields 9^ and 6'^, defined in terms of the 
auxiliary fields and ^ such that: the transverse part of 9^ is uniquely 
determined by the requirement that it should be the one that couples to the 
external current, and for the longitudinal part we choose it to be equal to 
the one of the field ^, i.e., 

0, = ^V^u^, - , (68) 

and 9'^ is instead defined in terms of $^ and ^ in such a way that the new 
action for the fields and 9'^ contains no mixing term. There are many 
consistent ways of achieve this, however giving the same bosonic action for 
the relevant field 9, but the simplest we found is 



0', = % - + ^V^IJ. (69) 



with 9^, of course, as given in (|68|). Note that, by virtue of (|69D, the longi- 
tudinal part of 9'^ is equal to the one of $, 

a-$ = d-9' (70) 

thus the functional delta function of $ becomes the gauge fixing for 6'. 

The relevant bosonic field will be 9^, since by virtue of the previous 
definitions 9'^ is decoupled and will consequently be ignored. The resulting 
generating functional for and the decoupled fermions is 
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+-e,[—{gC+-d^)V+ + —{gC^-d^)V-],.e, + iJ,e,,xdM ■ (71) 
We can also write it more explicitly, in the notation of ref. 0: 

Z{J) = Jvx^VxVe,expi^-Jd'x[xH^ + M)x + ^{d-9)' 



where 



F 



G 



(73) 



and F = F(— (9^), G = G{—d^) are given in terms of their Fourier transforms 
F and G by 



F(^2) 



m 



ATTk^ 



1 



( ) 2 



arcsin(l + 



fc2 



and 



G{k') = f + 



47r 27r I A: 



Aui^ 1 
arcsin(l + ^-75-)"^ 



(74) 



(75) 



where g is an arbitrary integer. 

Since x s^i^d x^ are decoupled from the relevant bosonic field 9, one recov- 
ers from (^) the bosonization of the Thirring model in 3 dimensions, with 
a bosonic action given, in the quadratic approximation of J($), by 

Sbos = J d'x[-^F^,{9)CiF,,{9) + '-9,C2e^Mx - ^-FlM + ^{d ■ 9f] 

(76) 

and with the correspondence between the currents 



^t(x)7^^(x) ^ ie^^xdu9x{x) . 



(77) 
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But the important new point is that, with our method, we are able to ex- 
press the initial fermion fields in terms of the bosonic ones and free fermionic 
fields, which was not the case in the usual approach to bosonization in 3 
dimensions, : 

^ = e-'^ {l - (^ + M)-i {^h^C^V^J, - ii^,{-d^)-'C2e^,xdJx+ n)}"' X 

(78) 

We end up this section by remarking that this redefinition of the fermionic 



fields yields the proper expectation value for the current (|77D , as can be 
verified by an application of (ESf) to this case. 
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4 The two ways of decoupling the fermions 
in 1 + 1 dimensions 



We shall here compare the usual way of decoupling massless fermions in 
d — 2, which relies upon the anomalous Jacobian method, and ours, which 
involves non-anomalous Jacobians, and moreover holds true in any number 
of dimensions (in particular d = 2). We perform the consistency check that 
both must lead to the same answer when applied to the two dimensional 
Abehan case. 

Let us begin by reviewing the anomalous Jacobian method. We start 
from the Euclidean generating functional 

Z{A) = jv^V^ eM- J d'^x^{^ + ie4)-^] . (79) 

The Dirac operator ^ + ie ^ may be written as 

^ + ie4 = ^ + ie ^{^^4) (80) 

where we understand as acting on A only, and not on the fermionic field 
that may eventually appear on the right. Similarly, ^ acts on (^~^ 4). 

Now, {^~^ /I) m 1 + 1 dimensions, leads in a natural way to the de- 
composition of the vector field A/^ into its longitudinal and transverse parts, 
since 

^-1/ = = a-2(a-^ + i75e^,a^A) (8i) 

or 

/ = - ^(^75^1 + ^2) (82) 

e 

where $1 = e d~'^e^j,dfj^Au and $2 = e d~'^d ■ A. 
Then one may write the Dirac operator as 

and achieve decoupling by performing the transformation 

^(x) = e-**2(^)-^^'^*i(^)x(a^) 

^(x) = ;^(a;)e**2(^)-^T5*i(^) . (84) 
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After the transformation ( ^4] ) is performed, the action for the new fermions 
is of course free 

SFi^,^;A) = 5f(x,X;0) (85) 

and there appears, due to the change of variables, an anomalous Jacobian 
J($i) 



J($i) = exp 



^ I d^x{d^<^i)'' = exp[|- I d?xF^,^F^J\ . (86) 
Thus, we have for Z{A) 

Z{A) = exp[|^ I SxF^^^F^,] I VxVx explj dPxx M ■ (87) 

Let us now show that our way of decoupling yields to an identical result. 
In Euclidean space, we must redefine the fields according to 

^(x) = {I +ie^'^ JiiY^^ x{x) 

m\x) = xHx){l + ie4^-Y^ , (88) 

after which the action becomes also free, as in the previous case. The differ- 
ence is that now we need to evaluate a (non anomalous) Jacobian 

_ det(^) 

•^^•^^ - ■ ^^^^ 



Again, due to the fact that we are considering the special case on two 
spacetime dimensions, we may use the important result that (^) receives 
a non-vanishing contribution only from the second order vacuum polarisa- 
tion, which therefore is the exact result for the logarithm of the Jacobian, 

JiA) = exp[^ I SxF^^^F^,] , (90) 

and we have 

Z{A) = exp[|^ j SxF^^^F^,] I Vx^Vx explj dPxx^ ^x] , (91) 
to be compared with 
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